The time-energy certainty relation 
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All energy measurements of a quantum system are prone to inaccuracies. In particular, if such 
measurements are carried over a finite period of time the accuracy of the result is affected by the 
length of that period. Here I show an upper bound on such inaccuracies and point out that the 
bound can be arbitrarily small if many copies of the system are available. 
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INTRODUCTION 



Since the onset of quantum mechanics the uncertainty 
principies have played the central role in attempts to 
comprehend this counterintuitive theory. The impor- 
tance of these principies springs from the search for lim- 
itations of the theory. They teli us what is possible 
and what is fundamentally unachievable in the quantum 
world. Although the mathematical structure of quantum 
mechanics is well understood, its physical interpretation 
and applicability to the phenomena of the observed world 
is still debated. In particular, there is some confusion 
about the status of various time-energy uncertainty rela- 
tions present in the textbooks on quantum theory P, Q . 
In this work I will focus on one important question that 
arises naturally in the context of the time-energy uncer- 
tainties: How accurate the energy measurement can be if 
it is carried over a small period of time? It has been no- 
ticed that there is no fundamental lower bound on such 
an accuracy no matter how short the measurement j3i0|· 
In the course of this Letter I derive an upper bound on 
the accuracy, which the energy of a quantum system can 
be determined with after a short measurement. To avoid 
further confusion on the subject, the derivation, which is 
simple enough, is presentcd in reasonable detail within 
the scope of this work. The final result is given in Eq. 



Lets consider a general setup where an evolution of 
a quantum system is governed by an unknown time- 
indepcndent Hamiltonian of discrete spectrum H. Any 
state of the system can be written in terms of the corre- 
sponding energy eigenstates |fc), but here I will consider 
only the states that are dominated only by a finite num- 
ber K of the nondegenerate eigenstates. Thus, the initial 



state of the system is expressed as 
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where the amplitudes ctfc, the eigenstates \k) and K are 
initially unknown. Later it will become clear that some 
states that are spanned over infinite number of the eigen- 
states can also be included into the consideration. 

I will be concerned only with the states at the discrete 
time steps t — nSt, n = 0, 1, N > K, that are denoted 
by \tp n ). In order to determine the eigenenergies hojk 
driving the evolution of \ip) one can measure the time 
autocorrelation function 
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with dk = | ctfe | 2 . This quantity can certainly be mea- 
sured since the wave function itself (not just the wave 
function modulus squared) can be measured with quan- 
tum tomography The eigenfrequencies are extracted 
from the eigenvalue problem for the evolution operator 
U, that evolves the system by one time step St, i.e. 

U\tpn) = l^n+l)) 



U\k) 



\k). 



(3) 



To fmd a matrix representation of the last equation it is 
convenient to arrange the first N states \ip n ) as columns 
of K x N matrix P. Inserting unity PP _1 before states 
\k) in Eq. ©, and multiplying by from the left results 
in the generalized eigenvalue problem 



u'm = 



-lUJ^St 



(4) 
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where U' = P^ÜP and S = P^P are N x N matri- 
ces with matrix elements expressed solely in terms of 
the autocorrelation function vàlues, U' mn = c„_ m +i and 
Smn — c n -m- The new eigenvectors are \k') = P _1 |fc). 

At this point the problem of extracting all the energy 
rclatcd information of the system in the initial state ji/'o) 



2 



is solved. Given the measured vàlues of the autocorrela- 
tion function c„, one finds the number K of eigenstates 
involved as a rank of the matrix S, eigenfrequencies from 
the solution of Eq. Q and coefficients dk by solving lin- 
ear set of equations with now known K and u>k ■ With 
known K, d k and uj k finding the expectation energy of the 
system or the corresponding Standard deviation is triv- 
ial. The method presented above goes under the name of 
harmònic inversion and its variations has been used by 
scientists on many occasions 0, 0, BEI • You can find an 
introduction to the method in Ref . [l0| . 



INACCURACIES 

The analysis, so far, did not account for the inevitable 
measurement inaccuracies. The experimental vàlues of 
the autocorrelation function differ from the unperturbed 
vàlues that enter Eq. 10) and this difference will affect the 
computed vàlues of the desired eigenfrequencies. The ef- 
fect of the experimental imperfectness will be further am- 
plificd by the requirement that the all measuring is done 
in a very short time interval T = NSt, which satisfies 
T < 2n/Auj max , with Au rnax = max^ fc {wj - u k } be- 
ing the largest possible frequency difference in Eq. J3J). 
In other words the total measurement time T is much 
shorter than the shortest characteristic time interval of 
the quantum evolution. 

Lets denote the measured value of the autocorrelation 
function by c„ = c n + r] n , where r/ n contains all the sta- 
tistical inaccuracies like noise and effect of finite number 
of the quantum system copies available. The systematic 
errors can be, in principle, eliminated and I assume they 
are absent. Also a time instant nSt, at which the mea- 
surement of c„ is performed, can be estimated with only 
finite precision, but in this work it assumes a precise value 
nSt, while the related inaccuracy is relegated to r) n . 

It is possible that even some states spanned over in- 
finite number of energy eigenstates can be considered. 
One would pick K < oo most important eigenstates i.e. 
with the greatest d^s and treat the remaining part as 
a noise contributing to rj n . However, this contribution 
would have a systematic character rather than statistical 
and should account for only a negligible part of the total 
error r) n . 

To summarize the discussion of sources of the inaccu- 
racies observe that all of them are contained in complex 
noise terms r\ n . In the absence of systematic errors, the 
noise level can be made arbitrarily small by increasing 
the number of copies of the quantum system undergoing 
the measurement. The time-energy uncertainty relation 
for a measurement on just one quantum system is intro- 
duced in ■ The largest magnitude of r) n is denoted by 

TJmax • 

In order to see the influence of the noise on the ex- 
tracted eigenfrequencies one must solve the generalized 



eigenvalue problem with corrupted by noise matrices 
U' and S. To do this lets first diagonalize the matrix S. 
S is a Hermitian non-negative defined N x N matrix. It 
has exactly K positive eigenvalues Ai > A 2 > ... > Xk > 
Xk+i = ... = Ajv = 0. The smallest positive eigenvalue 
will be denoted by X m in = Xk- The perturbed matrix S 
is still Hermitian and its eigenvalues differ from the un- 
perturbed ones by at most Nrj max , \X k ~ X k \ < N-q max . 
Thus, it is possible to teli the number of positive eigen- 
values of S only for limited level of noise 



< 



Xn 



2N 



(5) 



where the factor 1/2 comes from the fact that zero eigen- 
values are also perturbed up to Nrj max . This is the con- 
dition one must fulfill to extract all K freqüències that 
enter Eq. ©. 

In the next step one has to truncate the space spanned 
by eigenvectors of S only to those that correspond to K 
positive eigenvalues, so that the inverse of the truncated 
matrix S could be found. The eigenvectors to positive 
eigenvalues of S are arranged as a columns of N x K 
unitary matrix Q. The diagonal K x K matrix D = 
Q^SQ is nonsingular and can be inverted. Also the U 1 
is truncated to Q'U'Q. The Eq. Q takes the form an 
ordinary eigenvalue problem 



(6) 



In general, for an eigenvalue equation AX = XA, 
where X is a matrix of eigenvectors of A and A is the 
diagonal matrix of corresponding eigenvalues, the per- 
turbation SA of A results in perturbation 5A of A and 
X~ 1 (A+SA)X = A + SA. It follows that SA = X~ 1 5AX 
and using matrix 2-norm ||<5A|| < HA^H \\X\\ \\SA\\ = 
k(X)\\SA\\, with k(X) being the condition number of 
eigenvector matrix X [ 1 1| . 

Similarly, to estimate the influence of the perturbations 
on the eigenvalues in Eq. JüJ) lets rewrite the matrix on 
the left hand side as follows 



{D' 1 + SD-^Q^U' + SU')Q 
+5D- 1 &U'Q 



D- l &U'Q 

d- x qUu'q, (7) 



where SD^ 1 = Nij max dia,g(l/Xl, ...,1/X K ) and SU' has 
matrix elements rj n _ m+ i. The first term on the right 
hand side of Eq. J7J) has complex eigenvalues on the 
unit circle, which means that i-th row of QWQ must 
be of order of Ai to cancel the factor 1/Aj of D^ 1 . The 
second term must then have matrix elements at most of 
order of Nrj max /X m i n . The matrix elements of the third 
term are clearly of order of T) max /X m i n . The last two 
terms combined give an estimate of an error made when 
computing the eigenvalues of Eq. © 



-iu> k St i 



< k(P) 
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where k(P) is a condition number for matrix P and the 
fraction is an estimate of the Frobenius norm, the upper 
bound for 2-norm, of the two last matrices in Eq. J7J. 
The condition number for Q matrix is 1, because Q is 
unitary, and is skipped. 

The last inequality gives an estimate of frequency error 
made when evolution time step St is very short i.e. T <C 
27r/Aw maa . Namely, 



|íD fc - Lü k \5t < k(P) 



K(N + Ï) 

Am 



->)n 



(8) 



Multiplying both sides above by NTi one gets something 
opposite to a time-energy uncertainty relation, the upper 
bound rather than lower bound for the measured energy 
inaccuracies. This is why I refer to it as the time-energy 
certainty relation. The important remark is that the er- 
ror of the measured energies can be kept arbitrarily small, 
no matter how short the measurement was, by reducing 
the noise level rj max . It is assumed that the condition 
number k(P) is finite, as will be shown in the following. 

The certainty relation as presented in Eq. JSJl has littlc 
use, except perhaps a philosophical value, since it still 
contains rather cryptic quantities n(P) and X m i n . The 
prescription for calculating k(P) is provided in the next 
paragraph. 

Let as take a closer look at the K x N matrix P. As 
mentioned in the introduction the matrix columns are 
made of components of vectors \tp n ), Pkn = a k e~ tuJknSt . 
The procedure of finding unknown freqüències lük and 
corresponding amplitudes dk = \otk\ 2 starts from mea- 
sured vàlues of the time autocorrelation function (0). In 
the autocorrelation function the phases of the coefficients 
ak cancel out and only the modulai dk enter the Eq. 
Thus, in the definition of P one can safely use y/dk i n_ 
stead of ak and 



Pk 



d k e 



-iuJknSt 



(9) 



P is a rectangular matrix, in general, and its condition 
number is expressed in terms of singular vàlues rather 
then eigenvalues, k(P) = Pi/pk> where pi andp^ are the 
largest and the smallest singular vàlues of P, respectively. 
Notice that the matrix P has exactly K singular vàlues 
because there are exactly K linearly independent states 
\k) spanning the initial state in Eq. (JTJ. The singular 
vàlues pi are the square roots of the singular vàlues of 
Ptp = S. That is pi = y/%. Finally, 



k(P) 



Ai 



< 



ÍTt(S) 



A, 



(10) 



At that point the inequality JHJ becomes useful because 
the upper bound is expressed in terms of the trace and 
the smallest positive eigenvalue of the known matrix S. 
The eigenvalue is discussed in the following section. 



THE SMALLEST POSITIVE EIGENVALUE OF S 

The problem of the error estimates of the extracted 
freqüències is completed in Eq. JBJl and Eq. To 
better understand the relation between the possible error 
and the system parameters in Eq. lf2"jl. I will express the 
eigenvalue \ m in in terms of K, N , T, ujk and dk- 

It turns out that \ m in can be analytically estimated 
in the regime of short time span of the autocorrelation 
function. First helpful observation is that X m i n is a very 
small positive number in the limit where Aw max St -C 1. 
The S is Hermitian non-negative defined and it's matrix 
elements S mn — c„_ m approach ^DfcLi dk = T^(S)/N if 
the time step St is small. That means that the largest 
eigenvalue Ai approaches Tr(5) and there is little room 
left for other eigenvalues since the sum of all eigenvalues 
is Tr(S). Second hint comes from the property that the 
nonzero eigenvalues of the N x N matrix S = P^P are 
equal to nonzero eigenvalues of the K x K matrix PP> 
since both are squares of singular vàlues of the K x N 
matrix P. The last is a direct consequence of the singular 
value decomposition theorem fl3j |. 

Having established the smallness of \ m in and the 
equivalence of nonzero eigenvalues of P' P and PP\ one 
sees that the characteristic equation for S 



ajs[\ N + ... + ai\ + ao = 



fil) 



gives the estimate \ m in ~ — ff' wnere a o = det(PP^) 



and a% = — Yl,k=l ^ nor kk{PP^)- The matrix P has the 
following structure 



P = diag( v / dï, 



(12) 



where the n-th column of the KxN Vandermonde matrix 



V is given by the vector [ e -^i(n-i)St ) 
This leads to the determinant 



A" 



det(PP f ) = dct{VV r ) Y[ d k - 



-íüjk (n — l)5tl 



(13) 



fe=i 



The determinant of the product VV^ of the Vander- 
monde matrices in the limit of short time steps St is given 
by 



det(V-V-t) „ ( | 



A" 



A-l 



n 



l U 2 -r 



K-j 



K 



^) 2 



(14) 



The proof of this approximation is beyond the scope of 
this Lctter. It is easy, however, to numerically confirm 
correetness of Eq. ()14JI . 

A minor M.moikk(PP ) is also given by Eq. I|13|) where 
fc-th amplitude is dropped and fc-th frequency in Eq. i|14|) 
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is dropped and also K is replaced with K — 1 in both Eq. 
(O and Eq. 

Summarizing the calculations, the smallest eigenvalue 
of S is expressed in terms of unknown freqüències and 
amplitudes as follows 



(N + K-l)\ \{K-1)\ 
(N-K)\(2K-1) [(2K-2)\ 



(15) 



for K > 2, and 
for X = 2. 



(16) 



DISCUSSION OF THE RESULTS 

Despite the apparent complexity, the estimate of X m in 
has a simple structure. To see this, it is convenient to 
first define the effective frequency distance A 



number of freqüències K can blow the corresponding 
inaccuracies to a very high level. And this is perhaps 
the intuition that is hidden behind many time-energy 
uncertainty relations: "One cannot measure freqüències 
accurately in a very short time." This intuition does 
not reflect a fundamental principlc. It merely states the 
practical difhculty in analyzing systems involving many 
freqüències. 

Although there are time-energy uncertainty principies 
that are not based on an intuition but on a rigorous 
derivation from equations of quantum mechanics [ï^ . 
their meaning is different from the one considered in this 
work. 

In fact one could perform an experiment testing the 
absence of the lower bound on the energy measurement 
accuracy in a very short time, which is set by some time- 
energy uncertainty relations. Such an experiment would 
involve many copies of a two level quantum system pre- 
pared in the same state. Many copies are necessary to 
decrease the statistical noise rj max and the two energy 
levels would allow to keep X m in at a reasonable value. 
At fixed time interval T it is possible to adjust X m in, as 
long as X m i n ^ 1) by varying the number of time steps 
N, see Eq. |JT5J). 



di 



-J— = lv 



du 



■ UI i 
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for K > 2 and 

1 _ 1 

= 2did 2 (wi -w 2 ) 2 ' 

For dk = 1 the effective frequency distance 



K = 2. 



(17) 



(18) 



A e 



minjcij/j — u>j}, maxjwfe — ujj} 



with the tendeney to approach the left end of the inter- 
val, meaning that A characterizes the shortest frequency 
gaps. The dependence on the amplitudes dk suggests that 
A is dominated also by the weakest Fourier components 
of the autocorrelation function Eq. . 

The most important feature of X m i n is that X m i n oc 
(NStA) 2 ^- 1 ^. Using the total time span T — NSt one 
arrives at 



X mm oc (TA) 2(íf - 1} 



(19) 



Since only short measurements are considered 
TAu) max < 1 it follows that TA <C 1 and X m in 
decreases exponentially fast with the dimensionality of 
the Hilbert space the state lives in. That in turn 
implies that the upper bound on the measured frequency 
inaecuracy in Eq. JHJ is exponentially sensitive to the 
number of Fourier components K in Eq. Thus, for 
very short measurement times TA <C 1, even a moderate 



SUMMARY 

In this work the energy measurement on an evolving 
quantum system is considered. The evolution is governed 
by a discrete Hamiltonian and the initial state is spanned 
on a finite number K of significant energy eigenstates. 
The main result is the following upper bound 



líDfe - ui k \T < 



KN(N + l)yïr(S) 



X 



3/2 



(20) 



for inaecuracy \ü>k — u>k\ of measured frequency Qk with 
respect to the actual frequency u>k- The total measure- 
ment time T is short as compared to the characteristic 
time seales of the system, Taj ma x <§C 1, and is divided into 
N time intervals by the sampling measurements. If the 
noise r\ n is statistical then its maximal level rj max can 
be reduced by simultaneous use of many copies of the 
system. The systematic errors can, in principle, be elim- 
inated. Despite the short measurement time T the upper 
bound can be made arbitrarily low by the reduction of 
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